Itinerant Ferromagnetism for Mixed Valence Systems 
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We introduce a novel mechanism for the unusual itinerant ferromagnetism found in mixed valence 
systems like Ce(Rhi_a:Rui:)3B2 , Laa;Cei_2,Rh3B2 , US, USe, and UTe. With it we can provide an 
explanation for the long-unexplained large value of (~ 100° K) value and the maximum in the 
magnetization below Tc found experimentally. We also show that this novel itinerent ferromagnetism 
can be continuously connected with the localized case for which the energy scale is much smaller 
{Jrkky ^ 1°K). 
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Introduction. Understanding the mechanisms for fer- 
romagnetism usually involves a reconciliation of a local- 
ized electron picture, traced back to Heisenberg and 
an itinerant electron picture, traced back to Bloch [^. 
For many insulating materials, adding itinerant features, 
such as indirect exchange, to the localized picture brings 
satisfactory agreement with basic experimental features. 
For many metals, adding localized features, such as spin 
waves, to the itinerant picture has a similar effect. On 
the other hand it is now appreciated that for some novel 
classes of materials, such as the heavy Fermion and mixed 
valence materials, doctoring one picture or the other is 
often a questionable procedure. The electrons can display 
both localized and itinerant properties as these materials 
are neither good insulators nor good metals. 

In this letter we propose a simple mechanism for ferro- 
magnetism (FM) in heavy fermion mixed valence ma- 
terials. The mechanism relies on several specific en- 
ergy scales set by the band structure of these materi- 
als even though the electron correlation energy is the 
dominant energy scale. The simplicity and generality of 
the mechanism however should provide a useful frame- 
work for discussing ferromagnetism in broad classes of 
heavy fermion mixed valence materials often found in 
4/ and 5/ materials and their compounds. In fact we 
will show that it seems to explain some of the very un- 
usual long unexplained ferromagnetic properties of such 
Ce based compounds as Ce(Rhi_2.Ru^)3B2 Q and the 
uranian monochalcogenides B. 

The mechanism is based upon a two-band structure. 
The lower band near the zone center is dispersive and 
nearly free of electron correlation effects. Away from the 
zone center, this band is very flat and doubly occupied 
states experience strong correlation effects. The upper 
band is separated from the lower band by a gap that is 
smaller than the correlation energy. When the chemical 
potential is at the flat band level, the system becomes un- 
stable to the formation of a ferromagnetic state. There 
are compounds satisfying this mixed valence requirement 
under normal conditions or by either adding additional 
electrons (doping), applying pressure, or changing the 



temperature. This ferromagnetic state is a manifesta- 
tion of a Hund-like rule among electrons in band states 
as opposed to usual case of electrons in localized states 
(atomic orbitals). 

We will show that the periodic Anderson model (PAM) 
is a simple microscopic realization of this picture. Using 
quantum Monte Carlo simulations, we will demonstrate 
that the model admits a ferromagnetic ground state in 
the mixed valence regime with electron occupancies re- 
flecting our physical picture. We then develop a mean 
field picture that reproduces the ground state properties 
of the PAM and also allows us to compute finite tem- 
perature properties. Under certain conditions we find 
the finite temperature mean-field approximation predicts 
the same highly unusual behavior as a peak, below Tc, 
in the temperature dependence of the magnetization [||J|] 
and the deviation of the inverse susceptibility from Curie- 
Weiss behavior ||^. 

Model Hamiltonian. The PAM is described by the 
Hamiltonian 

^ = E + +Vj2idlf,, + 
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where d^- and create an electron with spin a \n d 
and / orbitals at lattice site i and = fl^fi^ is the 
number operator for the /-electrons of spin a at site i. 
The hopping amplitude t between d-orbitals is only to 
nearest-neighbor sites. The hopping amplitude V hy- 
bridizes different orbitals on the same site. 

When [/ = 0, the resulting non-interacting Hamilto- 
nian Hq has two dispersive bands: 



ek + e/± J(ek-e/)2-|-4U2 
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separated by an hybridization gap. For a chain, the band 
structure is illustrated in Fig. la. Irrespective of spatial 
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dimension, each band has regions where the k-states are 
dominantly of d or /-character. A mixed valence condi- 
tion arises when the Fermi energy Ep sits in a crossover 
region which is around the energy e j . Three energies are 
important: the hybridization gap A, the /-state disper- 
sion Sf = dE^ /dk\kp — hvp {vp is the Fermi velocity), 
and the Coulomb (correlation) energy U. (5/ is measur- 
ing the band dispersion close to the Fermi level. In the 
interacting problem we will assume that C/ > (5/, A. As 
we shall now discuss, the competition is between 5f and 
A which will be a competition between a ferromagnetic 
and paramagnetic state. 

Mechanism for Ferromagnetism. In Fig. la-c we illus- 
trate the physical mechanism for the ferromagnetic state 
which emerges when the electron filling p is arround 1 /4 
(one electron per lattice site) and the energy e/ of the 
/-orbitals is close to Ep. The latter condition defines 
the intermediate valence regime. 

To explain the mechanism, we first consider the non- 
interacting case (Fig. la). States in the lower band with 
mainly /-character have a small dispersion Sf because 
of the absence of direct hopping between the /-orbitals. 
In the interacting case, because the Coulomb interaction 
only affects the electrons in the /-orbitals, the electrons 
which were doubly occupying the states in the lower band 
with mainly d-character are practically unaffected. On 
the other hand, the electrons in the many states which are 
close to the Fermi level and have mainly an /-character 
will be strongly affected. These electrons spread to higher 
unoccupied k-states in the /-part of the band and polar- 
ize by an in-band Hund's rule (Fig. lb). To see the anal- 
ogy to Hund's rule more clearly, we consider the limiting 
case, represented in Fig. Ic, where Sf = and A ^ 0, re- 
ducing Fig. lb to a two level system in momentum space 
with each level being strongly degenerate (we only plot 
the states with predominant / character). If we add elec- 
trons to this two levels system it is easy to show that the 
ferromagnetic solution has the lowest energy and is there- 
fore the ground state. By polarizing, the spatial part of 
the wave function becomes antisymmetric and there is 
no double occupancy in the real space. In this way, the 
Coulomb repulsion is reduced to zero. The kinetic en- 
ergy has also the lowest possible value if the electrons 
are occupying the lower energy (degenerated) levels (see 
Fig. Ic). 

We now describe the conditions for the stability of the 
ferromagnetic state, and we will see how this mechanism 
departs from the one of Stoner |^ . To this end we state 
the condition U ^ A. From Fig. lb it is clear that the 
cost of the spreading of the ferromagnetic solution is an 
increase in the kinetic energy proportional to Sf. On the 
other hand, if we build up a nonmagnetic solution only 
using the states of the lower band, there will be a mini- 
mum delocalization for each particle due to the fact that 
a finite set of /-character states is in the upper band (see 
Fig. 1). This can be seen by constructing an / Wan- 
nier function with the lower band states. This Wannier 
function will have a delocalization lenght which depends 



on the k wave vector where the two bands are crossing. 
Therefore, by localizing in real space to avoid the double 
occupancy, the electrons will have a finite probability of 
occupying the /-states in the upper band. The energy 
cost per electron of occupancy is proportional to A. If 
the hybridization gap is much larger than the dispersion 
of the states, that is, A » Sf, the ferromagnetic state 
lies lower in energy than the nonmagnetic state. The 
energy of the excited nonmagnetic state is proportional 
to A. This can be easily seen from the limiting (and 
non-realistic) case Sf = ploted in Fig. Ic, for which 
the analogy with the atomic shells and the Hund's rule 
is evident. 
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FIG. 1. Illustration of the ferromagnetic mechanism. 

It is important to emphasize that the hybridization 
gap A and /-state dispersion Sf are the two basic in- 
gredients of our mechanism: A only appears if there 
is more than one band and Sf is small only if the hy- 
bridization is weak, li U > A, (5/, then A is the en- 
ergy scale of the paramagnetic state, while i5/ is the 
energy scale of the ferromagnetic solution. We can es- 
timate the magnitude of these scales for the PAM: for 
V < t/2, we get Sf ^ dek/dk\kpVy{e kp - e/)^ and 
A - V^/{2dt - e/) -I- v/(2di + e/)2 +41/2 _ 2dt. There- 
fore, if ef is close to the bottom of the lower band 
{ef + 2dt ^ t) and the Fermi level is close to the top 
{Ep ^ €f ~ V^/ {2dt — £/)), A is considerably larger than 
Sf and the ground state is ferromagnetic. We can see 
that the ferromagnetic solution is stable for comfortably 
realistic values of the parameters. It is clear from this 
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analysis the particle density p for a ferromagnetic solu- 
tion must be close to quarter filling. 

Quantum Monte Carlo Method. Using the constrained- 
path Monte Carlo method (CPMC), we computed the 
ground-state properties of the PAM on a square lattice. 
The CPMC method projects the ground state from an 
initial state \'iPt) by converting the iterative procedure 
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into a branched random walk. The details of the method 
are described elsewhere |7| as are our prior uses of it on 
the square PAM |,§. 

The defining characteristic of the CPMC method is its 
elimination of the fermion sign problem by excluding ran- 
dom walkers \(t>) that violate (?/'t|0) > 0. If [i/'t) were 
the exact ground state |'0o)) this constraint would gener- 
ate an exact elimination of the sign problem, and hence 
an exact solution. Extensive benchmarking indicates the 
CPMC method provides accurate estimates of the energy 
and various correlation functions. 

We prepared the trial state in specific values of the to- 
tal spin and z-component of spin 0] . Because the Hamil- 
tonian(^ conserves these quantum numbers, the itera- 
tive process produces a ground state with the same 
S and Sz- In the ferromagnetic regime, when the num- 
ber of lattice sites N was greater than the number of 
electrons iVg, the resulting energy as a function of S al- 
most always showed a minimum at a value of S satisfying 
< 5 < 5(2Af — Ne). Hence we typically found a par- 
tially polarized ferromagnetic ground state. 

For square lattices we also computed dependence of the 
static spin structure factor and the electron occupancies 
both on position and wave number space Most no- 
table for the present work is the wave vector dependence 
of the different spin components of the upper and lower 
band electron occupancies shown in Fig. 2. It has the 
same quantitative features as Fig. lb. 

Mean Field Theory. The quantum Monte Carlo simu- 
lations are limited to chains and square lattices of rela- 
tively small sizes at zero temperature. However, we have 
found that their predictions are well described by a sim- 
ple spin-polarized mean-field approximation. We can eas- 
ily extend this approximation to larger systems sizes in 
higher dimensions and at finite temperatures. We can use 
large lattices in three-dimensions at finite temperatures 
to compare with the results of experiments. 

The mean field Hamiltonian Hmf is 

Hmf = i?o + ^ + nl(n{^} - {nl){nL)) 

(4) 

(ri{| ) and (n^ ) are the determined in a self-consistent way 
within the mean-field theory. The translational symme- 
try of H implies that we can always find a translationally 
invariant ground state. This means that {n{^) and {n^^) 



will not depend on the site index i. Therefore there are 
two variational parameters to be determined from the 
mean-field equations: and 

Even though the mean field approximation can give 
a good description of the FM ground state, it is well 
known that such approximations overestimate the energy 
of the paramagnetic phase because they improperly esti- 
mate the real space correlations that are very important 
for the paramagnetic solution. Therefore it is crucial to 
check that the ground state is the FM one by using a 
more accurate method to evaluate the energy of the PM 
state. To this end we calculated the ground state of the 
PAM with 6x6 unit cells using the CPMC technique. We 
found excellent agreement with the energy and as seen in 
Fig. 2 the electron occupancies. 
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FIG. 2. Comparison between Mean Field and CPMC 
values of the occupation numbers of the non-interacting 
band states for each spin polarization in a 6 x 6 cluster 
(p = 7/24,1^ = 0.5t, U = t, and e/ = -it). On the a;-axis we 
only include the nonequivalent wave vectors ordered accord- 
ing to increasing non-interacting energies. 

Consequences. Bolstered by this agreement we used 
the mean-field approximation to calculate various ther- 
modynamic properties of the ferromagnetic phase for the 
three-dimensioanl cubic PAM. In Fig. 3a we show the 
temperature dependence of the magnetization for differ- 
ent values of ey. From the arguments given above it is 
clear that the magnetization comes from the electrons 
which are occupying the states with /-character. By in- 
creasing e f we are reducing the number of /-electrons and 
therefore the magnitude of the zero temperature magne- 
tization. Some critical value of e/ exists where the chem- 
ical potential starts departing from e / and a new energy 
scale £/ — Ep emerges in the system. This new energy 
scale is reflected in the appearance of a magnetization 
peak. (See the ef = —t case in Fig. 3a.) To understand 
the last statement one first has to realize that the zero 
temperature magnetization is small due to the reduction 
of the number of /-electrons. When the temperature is 
of the order of e/ — Ep, electrons are promoted electrons 
from the double occupied band states to the /-character 
states which have a large entropy (large density of states). 
The /-electrons will be polarized due to the energy con- 
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siderations above explained. In this way we can explain 
the origin of the magnetization peak below Tc- It is im- 
portant to note that the source for the large entropy is 
associated with charge and not with spin degrees of free- 
dom. This fact explains why a state with larger magne- 
tization has a higher entropy. The magnetization curves 
shown in Fig. 3a are in good qualitative agreement with 
the magnetization versus temperature data measured in 
Ce(Rhi_2:Ru2;)3B2 for different values of x H. 
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FIG. 3. Mean field results (p = 1/4,1/ = 0.5t, U = 2t) 
for (a) magnetization, (b) inverse susceptibility, (c) density of 
states and (d) specific heat. 

Another interesting aspect of this ferromagnetism is 
the deviation from linearity (Curie- Weiss behavior) for 
the inverse magnetic susceptibility above Tc over a large 
temperature range (see Fig. 3b) The behavior con- 
trasts that predicted by theories based on localized mag- 
netic moments. 

From above it is also clear that the change in entropy 
from the magnetic to the paramagnetic phase depends 
on the number for / electrons (spin degrees of freedom) . 
This dependence is seen in Fig. 3d where we plot the 
calculated specific heat for different values of e/. 

A characteristic of this ferromagnetic state that has 
consequences for the photoemission experiments is the 
implication that the quasi-particle dispersion should be 
close to that of the electrons in the non-interacting case. 
We recall that the ferromagnetic state described in Fig. 
1 is very similar to a non-interacting polarized state. As 
in the Stoner mechanism, the main role of the Coulomb 
repulsion is to polarize states with well defined momen- 
tum. The hybridization gap of the noninteracting solu- 
tion is most likely replaced by a pseudogap (see Fig. 3c). 
This pseudogap should be seen in the optical conductiv- 
ity measurements. 

Conclusions. Our quantum Monte Carlo simulations 
and mean-field calculations clearly show the existence of 
an itinerant ferromagnetic phase in mixed valence ma- 
terials that is supported by the PAM. This new phase 
describes qualitatively many experimental features such 
as an unusually large value for Tc(~ 100°K) and the max- 



imum in the magnetization below Tc PQ l that are found 
in such ferromagnetic compounds as Ce(Rhi_j:Ru2;)3B2 
§, Ce(Rhi_,Os,)3B2 i, and La.Cei-.RhgBa 0. 

Many previous theories of 4/ and 5/ electron materials 
treated these materials as systems of localized moments 
in the /-orbitals. These theories are thus unable to de- 
scribe an itinerant ferromagnetic phase and experimental 
consequences Q like the peak in the magnetization below 
Tc observed in some of these systems , the large value 
of Tc, the deviation of X~^{T) from the Curie- Weiss law 
above Tc in the uranium monochalcogenides H, and the 
mixed valence behavior of these compounds 0,|l^,^ . 

The physical picture just presented, when combined 
with our previous results allows a reconciliation of 
the localized and delocalized ferromagnetism pictures 
painted by Heisenberg and Bloch seventy years ago in 
the sense that it is possible to go continuously from the 
mixed valence FM state where the /-electrons are delo- 
calized to a FM state where there is one localized electron 
in each /-orbital |^ . In our picture one can do this by de- 
creasing e f from the Fermi level to a value near or below 
the bottom of the valence band. We note that the en- 
ergy scale of a localized ferromagnetic state is that of the 
RKKY interaction, which according to de Gennes's rule 
|lT| is of the order of 1°K, while the scale for the itin- 
erant case introduced here is that of the hybridization 
gap A, which is of the order of 100°K for heavy fermion 
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valence) to t 



In going from the delocalized (mixed 
le localized regime we thus expect a strong 
reduction of Tc accompanied by an increase of the zero 
temperature magnetization. (The T = magnetization 
is proportional to the / occupation number.) This ex- 
pectation is consistent with the observed behavior of the 
magnetization in Laa;Cei_:rRh3B2 as function of x pO[ |. 
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